A COMPUTER  SIMULATION  OF  A 
NAVXt  BOILER* 


xnx 


T.  Houlihan 


Approved  for  public  release;  distribution  unlimited 


UNCLASSIFIED 


-JWTV  CLASSIFICATION  OF  THIS  FAOI  (Wham  Dim  Kmtmrm4) 


REPORT  DOCUMENTATION  PAGE 

I • A* FONT  NUMBER 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


»•  RECIPIENT'S  CATALOG  NUNOCH 


1.  OOVT  ACCESSION  NO. 


«.  title  (mm  SuAmuj 

A Computer  Simulation  of  a Naval  Boiler 


S.  TYNE  OF  REPORT  S PERIOD  COVENEO 

Engineer's  Thesis; 
Septemhpr  1978 


t.  performino  ono.  repo* 

I.  CONTRACT  on  OAANT  NUMBER;.; 


T.  AUTHOR;  «J 

Wallace  Paul  Fini 


10.  PROORAM  ELEMENT.  PROJECT.  TASK 

ANSA  » WORK  UNIT  NUMBERS 


t.  PEAFOAMINQ  ORGANIZATION  NAME  AND  AOOAEU 

Naval  Postgraduate  School 
Monterey,  California  93940 


M.  CONTNOLLINO  OFFICE  NAME  AND  AOONESS 

Naval  Postgraduate  School 
Monterey,  California  93940 

14  MONI TORINO  AOENCV  name  A AOONCSVIf  lllmtmmt  I 


la.  REPORT  DATE 

September  1978 


IS.  NU« 


MLL 


i CmmUmlllm i OHt cm) 


IS.  SECURITY  CLASS,  (ml  Urn  r Span) 

Unclassified 


I la.  OECl  ASSlFICATION/OONNOAAOlNO 

schedule 


IS.  OISTNiauTION  STATEMENT  (ml  Urn  Rapart; 


Approved  for  public  release;  distribution  unlimited. 


IT.  OISTNiauTION  STATEMENT  (ml  Utm  aAalracf  mlmrmrn  Ik  DlmmM  70,  II  Ittmrmml  trmm  Rmpmn) 


IS.  SUPPLEMENTARY  NOTES 

IS.  KEY  WO  AOS  rCauiMi  an  ravaraa  m!4m  II  mmmmmmmrr  rnnrn  MmNO  Sr  Maa*  All  Sir; 

Navy  D-Type  Boiler 

Computer  Simulation 

Naval  Boiler 

Combustion  Engineering  Type  V-2M  Boiler 

10.  abstract  (Cmthm  mm  man  14m  II  nmmmmmmrr  mm4  lOmmtllr  Or  Olmmk  aalnj 

V*  A non-linear  computer  model  of  a U.S.  Navy 

D-Type 

Boiler  was  developed  using  lumped  parameters . 

The  program 

was  coded  in  IBM  CSMP-III  simulation  language. 

A companion 

routine  was  also  coded  to  permit  the  user  to  imDlement  the  —s. 

UNCLASSIFIED 

tuCUWTV  CL*MO'C»TIO«|  ow  Tmi 


(20.  ABSTRACT  Continued) 


■^“^model  for  a particular  boiler  using  only  readily  accessible 
data.  The  model  was  adapted  for  a Combustion  Engineering 
Type  V-2M  boiler  for  analysis. 


nr- 


□ □ k 


Approved  for  public  release;  distribution  unlimited 


A Computer  Simulation  of  a 
Naval  Boiler 

by 

Wallace  Paul  Fini 

Lieutenant  Commander,  United  States  Navy 
B.S.E.E.,  Worcester  Polytechnic  Institute,  1967 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degrees  of 


MASTER  OF  SCIENCE  IN  MECHANICAL  ENGINEERING 

and 

MECHANICAL  ENGINEER 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
September  1978 


3 :'y9  0 1 15  089 


ABSTRACT 


A non-linear  computer  model  of  a U.S.  Navy  D-Type 
boiler  was  developed  using  lumped  parameters.  The  program 
was  coded  in  IBM  CSMP-III  simulation  language.  A companion 
routine  was  also  coded  to  permit  the  user  to  implement  the 
model  for  a particular  boiler  using  only  readily  accessible 
data.  The  model  was  adapted  for  a Combustion  Engineering 
Type  V-2M  boiler  for  analysis. 
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NOMENCLATURE 


- Area  (ft  ) 

2 

Drum  Downcomer  Total  Cross  Sectional  Area  (ft  ) 

2 

Header  Downcomer  Total  Cross  Sectional  Area  (ft  ) 

2 

Main  Bank  Heat  Transfer  Surface  Area  (ft  ) 

2 

Main  Bank  Riser  Total  Cross  Sectional  Area  (ft  ) 

2 

Screen  Riser  Total  Cross  Sectional  Area  (ft  ) 
Becker's  Factor 
A General  Constant 

Specific  Heat  of  Liquid  (Btu/lbm/°F) 

Specific  Heat  of  Tube  Metal  (Btu/lbm/°F) 

Specific  Heat  (Btu/lbm/°F) 

Specific  Heat  of  Flue  Gas  in  the  Furnace 
(Btu/lbm/°F) 

Specific  Heat  of  Flue  Gas  at  the  Furnace 
Exit  (Btu/lbm/°F) 

Specific  Heat  of  Flue  Gas  at  Superheater  Exit 
(Btu/lbm/°F) 

Specific  Heat  of  Flue  Gas  at  Main  Bank  Exit 
(Btu/lbm/°F) 

- Constant  used  in  Rohsenow  Correlation 


C3-Cll 


Constants  Used  in  Various  Equations 
Diameter  (ft) 

Sum  of  the  Diameters  of  the  Drum  Downcomers 
(ft) 

Sum  of  the  Diameters  of  the  Screen  Downcomers 
(ft) 


9 


D - Sum  of  the  Diameters  of  the  Main  Bank  Tubes 

r (ft) 

D - Sum  of  the  Diameters  of  the  Screen  Tubes  (ft) 

s 

EB  - Total  Heat  Input  to  the  Furnace  (Btu/sec) 

EDI  - Heat  Transferred  from  the  Desuperheated 

Steam  to  the  Tubes  (Btu/sec) 

ED2  - Heat  Transferred  from  the  Superheater  Tubes 

to  the  Steam  (Btu/sec) 

ER1  - Heat  Transfer  Rate  from  the  Flue  Gas  to  the 

Screen  Tubes  (Btu/sec) 

ER2  - Heat  Transfer  Rate  from  the  Flue  Gas  to  the 

Superheater  (Btu/sec) 

ER3  - Heat  Transfer  Rate  from  the  Flue  Gas  to  the 

Main  Bank  (Btu/sec) 

ER4  - Heat  Transfer  Rate  from  the  Flue  Gas  to  the 

Economizer  (Btu/sec) 

E - Heat  Transfer  Rate  to  the  Riser  Liquid 

wr  (Btu/sec) 

EWl  - Heat  Transfer  Rate  from  the  Economizer  Tubes 

to  the  Feed  Water  (Btu/sec) 

EW2  - Heat  Transfer  Rate  from  the  Desuperheater 

Tubes  to  the  Drum  (Btu/sec) 


EW3 

EW4 


Heat  Transfer  Rate  to  the  Screen  Riser  Liquid 
(Btu/sec) 

Heat  Transfer  Rate  to  the  Main  Bank  Riser 
Liquid  (Btu/sec) 

Friction  Factor 

Drum  Downcomer  Friction  Factor 
Heater  Downcomer  Friction  Factor 
Main  Bank  Riser  Friction  Factor 
Screen  Riser  Friction  Factor 

2 

Acceleration  of  Gravity  - 32.2  ft/sec 
Mass  Flow  Rate  (lbm/sec) 
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GB 

GDI 

GD2 

GD3 

GD4 

GL 

GR 

GS 

Gr 

Gw 

Gjri 

Gwr 

Gws 

HDl 

Hfg 

HO 

HNB 

HNBR 

HNBS 

HWl 

HW2 

HW3 

HW4 

HW6 


Fuel  Flow  Rate  (lbm/sec) 

Total  Steam  Flow  Rate  fron  Boiler  (lbm/sec) 

Desuperheated  Steam  Flow  Rate  (lbm/sec) 

Steam  Condensing  from  Vapor  to  Liquid  in  the 
Steam  Drum  (lbm/sec) 

Superheated  Steam  Flow  (lbm/sec) 

Air  Mass  Flow  Rate  (lbm/sec) 

Total  Flue  Gas  Flow  Rate  (lbm/sec) 

Screen  Riser  Circulation  Flow  Rate  (lbm/sec) 

Riser  Circulation  Flow  Rate  (lbm/sec) 

Downcomer  Circulation  Flow  Rate  (lbm/sec) 

Main  Bank  Riser  Mass  Flow  Rate  (lbm/sec) 

Drum  Downcomer  Circulation  Flow  Rate 
( lbm/sec) 

Header  Downcomer  Circulation  Flow  Rate 
( lbm/sec) 

Enthalpy  of  Saturated  Steam  (Btu/lbm) 

Heat  of  Vaporization  (Btu/lbm) 

Heating  Value  of  Fuel  Oil  (Btu/lbm) 
Non-boiling  Height  (ft) 

Main  Bank  Riser  Non-boiling  Height  (ft) 

Screen  Riser  Non-boiling  height  (ft) 

Enthalpy  of  Saturated  Water  (Btu/lbm) 

Enthalpy  of  Downcomer  Liquid  (Btu/lbm) 

Enthalpy  of  Feed  at  Economizer  Inlet 
(Btu/lbm) 

Enthalpy  of  Feed  at  Economizer  Outlet 
(Btu/lbm) 

Enthalpy  of  Water  in  Lower  Drum  (Btu/lbm) 
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Hxr 

Hxs 

k 

K 

Ks 

KR1-KR4 

l 

L 

LB 

Ldr 

Lds 

Ldlr 

Ldls 

Lr 

Lrl 

Ls 

Lsl 

LMTD 

M 

MMl 

MM2 

MM3 

MM4 

MM5 

MR1 

MW1 


Effective  Height  of  Main  Bank/Drum 
Downcomer  Circuit  (ft) 

Effective  Height  of  Screen/Header  Downcomer 
Circuit  (ft) 

Thermal  Conductivity  (Btu/hr-f t-°F) 

a general  constant 

Equivalent  Sand  Roughness  (in) 

Constants  Used  in  the  Gas  Heat  Transfer  Equations 
Length  (ft) 

Tube  Length  (ft) 

Boiling  Length  (ft) 

Total  Length  of  Drum  Downcomer  Tubes  (ft) 

Total  Length  of  Header  Downcomer  Tubes  (ft) 

Avg.  Length  of  Drum  Downcomer  Tubes  (ft) 

Avg.  Length  of  Header  Downcomer  Tubes  (ft) 

Total  Length  of  Main  Bank  Tubes  (ft) 

Avg.  Length  of  Main  Bank  Tubes  ‘(ft) 

Total  Length  of  Screen  Tubes  (ft-) 

Avg.  Length  of  Screen  Tubes  (ft) 

Log  Mean  Temperaturfe  Difference  (°F) 

Mass  (lbm) 

Mass  of  Screen  Tube  Metal  (lbm) 

Mass  of  Superheater  Tube  Metal  (lbm) 

Mass  of  Riser  Tube  Metal  (lbm) 

Mass  of  Economizer  Tube  Metal  (lbm) 

Mass  of  Desuperheater  Tube  Metal  (lbm) 

Mass  of  Flue  Gas  in  Furnace  (lbm) 

Mass  of  Water  in  Steam  Drum  (lbm) 


p 


Pressure  (psia) 

Pressure  of  Steam  Volume  in  Steam  Drum  (psi) 
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PD 

PD1 

PD2 

PD2 

Pr 

Psat 

PW1 

PW2 

q 

R 

Ral 

Ra2 

RD1 

Re 

RW 

RS 

Rtf 

t 

T 

Tsat 

TD 

TDl 

TD2 


Saturation  Pressure  Corresponding  to  Enthalpy 
of  Water  in  the  Steam  Drum  (psi) 

Pressure  at  Superheater  Outlet  (psia) 

Pressure  at  Desuperheater  Outlet  (psia) 

Prandtl  Number 

Saturation  Pressure  (psia) 

Pressure  in  Water  Drum  (psi) 

Pressure  in  Screen  Header  (psi) 

Heat  Transfer  Rate  (Btu/sec) 

Tube  Radius  (in) 

Average  Density  in  the  Main  Bank  Risers 
(lbm/ft3) 

- Average  Density  in  the  Screen  Risers 
(lbm/ft3) 

3 

Density  of  Saturated  Steam  ('lbm/ft  ) 

Reynolds  Number 

Density  of  Saturated  Water  (lbm/ft3) 

Density  of  Two-Phase  Mixture  at  Screen 
Outlet  (lbm/ft3) 

- Density  of  Two^Phase  Mixture  at  Main  Bank 
Outlet  (lbm/ft  j 

Time  (sec) 

Temperature  (°F) 

Saturation  Temperature  (°F) 

Temperature  of  Steam  Drum  Vapor  (°F) 

Temperature  of  Saturated  Steam  Corresponding 
to  Drum  Enthalpy  (°F) 

- Superheater  Outlet  Temperature  (°F) 
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TD3 


I 


! 


I 


TM1 

TM2 

TM3 

TM4 

TM5 

TR1 

TR2 

TR3 

TR4 

TR5 

TRa2 

TRa3 

TRa4 

TW3 

TW4 

TO 

Ub 

Um 

V 

VB 

VD1 

VENT 

VNDR 

VR 

VS 


Desuperheater  Outlet  Temperature  (°F) 

Screen  Metal  Temperature  (°F) 

Superheater  Metal  Temperature  (°F) 

Main  Bank  Metal  Temperature  (°F) 

Economizer  Metal  Temperature  (°F) 

Desuperheater  Metal  Temperature  (°F) 

Temperature  of  Flue  Gas  in  Furnace  (°F) 

Temperature  of  Flue  Gas  Leaving  Furnace  (°F) 

Temperature  of  Flue  Gas  Leaving  Superheater  (°F) 

Temperature  of  Flue  Gas  Leaving  Main  Bank  (°F) 

Temperature  of  Flue  Gas  Leaving  Economizer  (°F) 

Average  Flue  Gas  Temperature  Across  the 
Superheater  (°F) 

Average  Flue  Gas  Temperature  Across  the 
Main  Bank  (°F) 

Average  Flue  Gas  Temperature  Across  the 
Economizer  (°F) 

Feed  Inlet  Temperature  (°F) 

Feed  Outlet  Temperature  (°F) 

Ambient  Temperature 

Bubble  Rise  Velocity  (ft/sec) 

Mean  Velocity  (ft/sec) 

Volume  (ft3) 

Boiling  Volume  (ft3) 

3 

Volume  of  Steam  in  Steam  Drum  (ft  ) 

Percent  Throttle  Valve  Opening 
Water  Drum  Volume  (ft3) 

Main  Riser  Bank  Fluid  Volume  (ft3) 

Screen  Riser  Bank  Fluid  Volume  (ft3) 
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J 


V 


X 

X0 

Xs 

Xe 

Z1 

Z2 

Z 3 

Z4 

Z5 

Z6 

Z7 

Z8 


3 

Specific  Volume  (ft  /lbm) 

Quality  (percent) 

Quality  at  Main  Bank  Riser  Outlet  (percent) 
Quality  at  Screen  Riser  Outlet  (percent) 
Quality  at  Riser  Exit  (percent) 

Drum  Downcomer  Bend  Loss  Factor 
Drum  Downcomer  Exit  Loss  Factor 

- Main  Bank  Riser  Bend  Loss  Factor 
Main  Bank  Riser  Exit  Loss  Factor 
Screen  Riser  Bend  Loss  Factor 
Screen  Riser  Exit  Loss  Factor 

- Header  Downcomer  Bend  Loss  Factor 
Header  Downcomer  Exit  Loss  Factor 


Greek  Letters 


u 

p 

o 


pDal 

pDa2 


- Viscosity  (lbm/ft-sec) 

Density  (lbm/ft3) 

- Surface  Tension  (lbf/ft) 

- Average  Density  in  Superheater  (lbm/ft3) 

- Average  Density  in  Desuperheater  (lbm/ft3) 


Subscripts 


a 

- Average 

Liquid 

g 

- Vapor 

M,m 

- Metal 

r 

s 

dr 

ds 

D 

R 

W 


Far  Field  Average 

- Main  Bank 
Screen 

- Drum  Downcomer 
Header  Downcomer 
Steam 

Flue  Gas 

- Water 
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I . INTRODUCTION 


The  increasing  attractiveness  and  efficiency  of  computer- 
aided  design  techniques  have  established  a need  for 
reasonably  accurate  computer  models  of  all  types  of 
engineering  systems.  The  design  of  multivariable  control- 
lers, in  particular,  requires  the  availability  of  state 
information  about  a system  which  is  not  available  from 
conventional  transfer  function  simulations.  A number  of 
computer  models  of  marine  boilers  have  been  proposed 
[1,2, 3,4].  All  of  these  efforts  share  one  or  more  of  the 
following  features: 

a.  Lumped  parameters  are  used. 

b.  Bubble  formation  effects  with  respect  to  shrink  and 
swell  are  not  modeled. 

c.  The  furnace  screen  and  main  riser  bank  are  lumped 
together  as  are  their  respective  downcomers  [1,3,4]. 

d.  The  model  is  linearized  using  a Taylor  expansion 
or  perturbation  technique  [3] . 

Use  of  lumped  parameters  is  dictated  by  the  need  to  maintain 
computational  efficiency,  particularly  when  the  model  is 
used  in  control  system  design  and  optimization  schemes. 

In  the  past,  effects  due  to  vapor  formation  in  the  drum 
have  not  been  included  due  to  the  lack  of  a suitable  model. 

A simple  model  based  on  the  Rohsenow  correlation  [5]  was 
developed  for  use  in  the  current  simulation.  However,  its 
contribution  to  shrink  and  swell  effects  was  negligible  and 
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it  was  eliminated  from  the  equation  set.  Regarding  the 
tube  banks,  Haraguchi  and  Akasara  [2]  split  out  the  various 
riser  and  downcomer  loops,  including  intermediate  headers. 
Considering  the  general  level  of  uncertainty  inherent  in 
the  correlations  used  in  this  model,  this  appears  to  be  an 
excessive  degree  of  complication.  A compromise  position 
was  taken  in  the  present  effort  in  that  the  screen  and  main 
riser  bank  and  their  associated  downcomers  were  modeled 
separately,  but  intermediate  headers  within  the  screen 
circuit  were  not  addressed  individually  since  their  configura- 
tion varies  widely  from  one  boiler  to  another. 

The  objective  of  this  thesis  was  to  develop  a general, 
non-linear  model  of  a Navy  boiler  for  use  in  simulations 
covering  the  full  operating  range  of  the  propulsion  plant. 
However,  since  a linearized  version  is  normally  needed  for 
multivariable  controller  design,  CENTRL,  a computer  code 
developed  locally  for  use  with  CSMP  models,  could  be  used 
to  generate  linearized  state  space  representations  described 
by  the  matrix  equations: 

Si}  * L»l  SA*1  '-[M 

Lc  1 


This  approach  would  produce  linearized  representations 
covering  the  entire  boiler  operating  range  in  a piecewise 
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fashion,  rather  than  attempting  to  extend  the  range  of 
applicability  of  a single  state  space  or  linear  model. 

Since  this  author  found  the  notation  used  in  Reference 
[1]  particularly  readable,  it  was  adopted  for  use  in  this 


document . 


II.  DESIGN  CONSIDERATIONS 


A.  D-TYPE  BOILER 

The  single  furnace  (D-Type) , uncontrolled-superheat 
boiler  was  chosen  as  the  basis  for  the  model  since  it  is 
the  predominant  design  in  the  current  U.S.  Navy  inventory. 
The  nominal  operating  pressures  range  from  600  to  1200  psi, 
and  all  designs  are  conceptually  similar.  A basic  schematic 
is  shown  in  Figure  1.  These  boilers  are  typified  by 
relatively  fast  response  times  for  use  in  high  speed 
maneuvering  vessels.  This  puts  heavy  emphasis  on  properly 
designed  control  systems.  Manual  control  of  these  boilers 
is  essentially  infeasible  except  under  the  most  favorable 
circumstances,  e.g.  steady  steaming  at  low  speeds.  Current 
conventional  control  systems  used  by  the  Navy  are  of  the 
same  generic  type  with  relatively  minor  customization  except 
in  hardware  from  one  installation  to  the  next.  These 
designs  have  been  based  on  analog  models  intended  to  provide 
reasonable  dynamic  responses  for  pressure,  water  level,  and 
steam  flow  which  are  the  only  state  parameters  measured  by 
conventional  controllers.  Although  these  control  systems 
have  in  general  proven  satisfactory,  there  is  some  indica- 
tion [6]  that  multivariable  control  systems  can  effect 
significant  improvement  in  both  response  time  and  demands 
on  auxiliary  equipment. 


A typical  boiler  of  this  type,  the  Combustion  Engineering 
Type  V-2M  is  shown  in  Figure  2.  The  main  components  of  the 
boiler  are  the  furnace  which  is  surrounded  on  all  sides 
except  the  floor  by  screen  tubes,  the  superheater,  main 
riser  banks,  economizer,  water  drum,  steam  drum,  desuper- 
heater, downcomers,  and  screen  headers.  Navy  boilers 
do  not  include  air  heaters  and  the  only  air  preheating 
which  occurs  is  accomplished  by  passing  the  incoming  air 
over  the  exhaust  stack  liner.  This  procedure  generally 
produces  an  air  temperature  of  about  100-130°F  for  an 
ambient  of  70°F.  In  contrast  to  the  boilers  used  as  the 
basis  for  previous  models,  the  desuperheater  is  located 
in  the  water  drum  rather  than  in  the  steam  drum.  All  of 
the  steam  generated  passes  through  the  superheater  prior 
to  being  split  out  as  main  superheated  steam  and  auxiliary 
desuperheated  steam. 

All  of  these  boilers  are  of  the  natural  circulation 
type.  The  "throttle  valve"  was  considered  to  be  in  the 
superheated  main  steam  line. 

B.  MODEL  FLEXIBILITY 

The  model  was  designed  to  be  as  flexible  as  possible 
to  permit  its  relatively  simple  adaptation  to  a particular 
boiler.  Since  the  boiler  technical  manual  or  equivalent, 

I 

the  steam  and  gas  tables,  and  standard  engineering  hand- 
books may  be  presumed  to  be  available,  the  model  was  designed 
to  utilize  data  only  from  these  sources  for  implementation. 
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To  convert  this  input  data  into  the  form  required  by  the 
model,  an  interface  program  was  developed  to  adjust 
parameters  and  determine  initial  conditions  at  the  operating 
point  for  which  the  input  data  were  extracted.  This 
approach  was  taken  to  avoid  the  necessity  of  redesigning 
the  model  for  each  operating  range  or  specific  boiler 
design  being  considered. 

C.  ASSUMPTIONS 

To  permit  tractable  computation,  a number  of  simplifying 
assumptions  were  made  in  the  formulation  of  the  model; 

1.  Feed  water  entering  the  drum  mixes  thoroughly  with 
the  drum  liquid. 

) 

2.  No  steam  generation  occurs  in  the  steam  drum  except 
through  mass  transfer  at  the  water-vapor  interface. 

3.  There  is  a uniform  distribution  of  heat  flux  in  the 
generating  tubes. 

4.  The  mixture  of  steam  and  water  in  the  generating 
tubes  is  homogeneous. 

5.  Quality  varies  linearly  in  the  generating  tubes. 

6.  Heat  losses  to  the  ambient  except  through  the 
exhaust  gases  are  neglected. 

7.  All  downcomers  and  risers  in  a particular  circuit 
have  the  same  length  as  the  average  length,  and  the 
shapes  and  diameters  are  equal.  Also  the  number  of 
tubes  is  the  same  as  in  the  actual  boiler. 

8.  The  slip  ratio  is  1.0  (no  slip). 
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Since  the  feed  enters  through  a long,  horizontal  feed 
pipe,  and  the  drum  liquid  is  in  a state  of  turbulence, 
the  first  assumption  is  well  justified.  Although  the 
second  assumption  is  not  strictly  correct  during  transients 
when  there  is  a significant  change  in  drum  pressure,  and 
hence  saturation  temperature,  the  effect  of  steam  generated 
at  the  drum  wall-liquid  interface  appears  to  be  negligible 
based  on  an  analysis  of  this  factor  in  connection  with 

this  study.  The  assumption  of  uniform  heat  flux  is  ' 

1 

recommended  by  Reference  [7]  as  a reasonable  alternative 
to  a stepwise  integration  over  the  tube  height.  The 
assumption  of  a linear  change  in  quality  is  a consequence 

of  the  assumption  of  a linear  heat  load.  For  a well  y 

t 

designed  boiler,  heat  losses  through  the  casing  are  small 
compared  to  the  total  energy  generated.  The  effects  of 
any  such  losses  will  be  reflected  in  small  errors  in  the 
computed  heat  capacities  of  the  gas  and  in  the  heat  transfer 

coefficients.  The  assumptions  regarding  tube  geometry  are  i 

made  primarily  for  computational  simplification.  The 

last  assumption  regarding  slip  is  not  an  essential  one, 

and  higher  slip  ratios  could  easily  be  considered  if 

necessary  to  the  analysis  of  a particular  boiler  design. 

In  addition  to  these  explicitly  stated  assumptions,  those 
inherent  in  the  correlations  employed  in  the  model  formula- 
tion are  also  considered  to  have  been  made. 

■it 
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III.  MODEL  ANALYSIS 


A.  GENERAL 

The  model  was  envisioned  as  an  assemblage  of  control 
volumes  for  which  the  rates  of  energy  transfer  are 
described  by  appropriate  equations  to  obtain  the  rate  of 
change  of  energy  storage.  The  resulting  differential 
equations  were  integrated  to  obtain  the  "state"  of  the 
system.  In  this  way,  the  states  can  be  chosen  to  have 
specific  physical  interpretations.  The  system  is  shown 
diagrammatically  in  Figure  3. 

The  following  control  volumes  were  considered  in  the 
model  formulation: 


1. 

Riser 

tube  walls 

2. 

Superheater  tube  walls 

3. 

Economizer  mass 

4. 

Desuperheater  tube  walls 

5. 

Riser 

fluid  (kinetic  energy) 

6. 

Riser 

fluid  (thermal  energy) 

7. 

Steam 

drum  water  volume 

8. 

Steam 

drum  steam  volume 

9. 

Water 

drum 

10. 

Furnace 

B.  FLUE  GAS  AND  COMBUSTION  SYSTEM 

The  flue  gas  and  combustion  system  is  shown  schematically 
in  Figure  4.  Heat  is  carried  into  the  furnace  through  the 
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sensible  heat  of  the  air  and  fuel,  and  the  heat  released 
in  combustion.  Since  the  fuel-air  ratio  is  relatively 
constant,  and  the  fuel  and  air  have  relatively  constant 
inlet  temperatures,  the  sensible  heats  were  lumped  with 
the  heating  value  of  the  fuel.  The  mass  flow  rate  into 
the  furnace  is  given  by 

G-r.  * +•  G\_ 

and  the  total  heat  input  is 

^•a  3 Gii  * Ho  . 


The  heat  transfer  to  the  furnace  screen  is  primarily  through 
radiation,  i.e. , 


lit 


joo 


'fai  * 4io 
/ oo 


The  heat  balance  for  the  mass  of  flue  gas  in  the  furnace  is 


CfciTat ) = Es  - £*4  - £* 


where  the  last  term  represents  the  heat  carried  off  by  the 
flue  gas.  Heat  transfer  to  the  remaining  tube  banks  is 
primarily  through  convection.  The  heat  given  up  by  the 
flue  gas  is  given  by  the  general  relation 


and  the  heat  transfer  to  the  tubes  may  be  estimated  with 
the  Grimson  Correlation  [8] 


. 


I 


Pr  ^ ~ "Gi ) 

Since  the  physical  properties  in  the  above  equation  do 
not  change  significantly  over  the  expected  range  of 
temperature  variations,  the  correlation  is  used  in  the 
form 

Tw  » K GfT  C - Try,  ) 

Since  the  tube  geometry  is  not  known  in  advance,  an  n = 0.6 
was  used  as  a reasonable  average  for  most  geometries  found 
in  practice.  The  equations  for  the  tube  banks  were: 


ER.i.  « Gg.  Caia.  (Ta.z  -"Tm) 

Ea.2.  * Ga.'u  -T»ax,) 

To.*!*  v To-0  /to 


2 . Main  Bank 

Ew  m G«.Cdz.  (Taj  — Tm) 

* K*s6*'U(Ta*j  -Tws) 

Ctr.»  «►  T«.m)  1 2-o 


3 . Economizer 

- Gr.  Cfci  (Te^  -Te.s ) 

CT(i&.^  — Timh) 

Ta-M  « + Tr.*')  / a.o 

i ' 

The  arithmetic  average  temperature  for  the  gas  was  used  in 

. 

the  driving  potential  AT  rather  than  the  log  mean 
temperature  differences  ( LMTD ) , because  a comparison  of 
the  results  obtained  for  the  two  methods  indicated  the 
differences  were  slight  for  the  temperatures  involved  and 
the  simpler  average  was  more  computationally  efficient. 

C.  WATER-STEAM  CIRCULATION  SYSTEM 

As  is  shown  in  Figure  5,  the  basic  circulation  system 
consists  of  the  steam  drum,  water  drum  or  header,  and  the 

. 

riser  and  downcomer  tubes.  The  driving  force  for  circulation 

. 
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is  the  difference  in  density  between  the  riser  and  down- 
comer loops.  The  pressure  loss  around  the  circuit  balances 
this  driving  force.  For  the  downcomers,  the  force  balance  is 


for  the  drum  downcomers , and 


for  the  screen  downcomer. 

For  the  risers,  the  equations  are: 


Pwi  - +z,J£L+liJ£ 1_ 

+ e.l3«»*+^^,T  (7?) 

for  the  screen  and  main  bank,  respectively.  The  friction 
factors  are  calculated  from  [9] 


(/■  7y  -2lojg/n%)  J 
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The  friction  factors  for  the  risers  are  multiplied  by 
Becker's  factor  [10]  to  account  for  two-phase  effects 


Be  - i *■  /uoyg.sg(^-^ 


9i 


The  relationship  between  riser  and  downcomer  mass  flow 
rates  is  given  by  the  continuity  equation 


_L  f 

it 


Gw"  Gy 


The  average  density  in  the  risers  is  found  from 


--f 


From  the  assumption  that  the  quality  varies  linearly,  the 
density  may  be  obtained  as  a function  of  Z from 


!(»= 


V(  * L 


Regrouping  the  terms  we  get 


f(i)  - 
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This  expression  can  be  integrated  in  closed  form  to  yield 


(ftOJi 


J. 


-k-  In 


^ v*+  v* 


Combining  terms  and  simplifying  finally  yields 


( fi«IJ " xh; u ■ 

V*  4 L 


.fo  Vfr 


Vi 


+ 1 


Prom  this  result,  p is  obtained  as 

cl  V 


1 

La 

W 

XcVw  tl 

L 

T 

A 

# 


An  energy  balance  on  the  riser  tube  banks  yields  the 
quality 

it  { f„  ( ««  * -£■  V, 

s Wwii  “ Grt  ^Hwi  f V E*y\r  Grot  (jtwi 
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where  VB  is  the  boiling  volume,  and  the  last  term  on  the 
right  accounts  for  heat  added  in  the  non-boiling  height 
to  raise  the  enthalpy  to  saturated  conditions. 

The  non-boiling  height  is  obtained  from 

u a*.  L — 

where  qg  is  the  sensible  heat  added  to  the  liquid  in  the 
non-boiling  height,  and  qfc  is  the  total  heat  input.  The 
ratio  q /q.  is  given  by 

Hwi**  Hwr 

(vU.  8*3)-  H-* 

The  boiling  volume  is  related  to  the  non-boiling  height 
by 

vB  = L~  ^ V ■ 

u 


With  substitution  of  the  appropriate  parameters,  the  above 
equations  were  applied  to  both  circulation  loops. 

The  mass  balance  of  the  water  in  the  steam  drum  is 


it 


* Gro  ^ G 6 S ► & VX  ~ G tvy  — Gn*>S 
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and  the  energy  balance  is  given  by 

Mwl)  m Go  (l-X.)  rtwt  Gi  C '“■*»)  Hwi 

+■  Goj  H +■  Gwi  ~ Gu»r  tfm.  - Gu>*  |4m»x 

Since  all  the  incoming  feed  must  be  passed  through  the 
riser  tubes  to  be  heated,  an  energy  balance  on  the 
downcomers  gives 

Gvox  ^ f &wr  4*  Gwi  — Quit. 

J 

* C 6w»"  t HwX 

' V 

Gwi  (,Wu»^  — rtwt ) v +■  0*>&  ) 14  wi 

s C ^***r  +■  Hu>i 

| 

Solving  for  HW2,  the  downcomer  enthalpy,  we  obtain 

IW  * S^TS,  -K“*)  v «“*■ 
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D.  WATER-STEAM  HEAT  TRANSFER 


1.  Risers 

The  heat  transferred  to  the  riser  tubers  is  used  to 
bring  the  incoming  flow  to  saturation  conditions  in  the 
non-boiling  height,  and  then  generate  steam  in  the  remain- 
ing portion  of  the  tube.  To  permit  a tractable  equation 
set,  however,  the  tube  walls  were  assumed  to  be  at  a uniform 
temperature,  and  all  the  heat  was  assumed  to  have  been 
transferred  under  fully  developed  flow  boiling  conditions. 
Since  the  non-boiling  height  for  the  main  bank  was  found  to 
be  small,  the  assumption  was  reasonable  for  that  circuit. 

The  non-boiling  heights  were  higher  in  the  screen  circuit; 
however,  the  error  in  temperature  introduced  was  less  than 
three  percent  which  was  deemed  an  acceptable  level.  For 
fully  developed  flow  boiling,  Tong  [11]  recommends 


. _ fi. 


i.m  * 10* 


Therefore  the  heat  transfer  to  the  risers  was  given  by 


Ew3  - Kwj  Pox/a  (-r*i  -Toi  )J 

- K-M  Pol'</,  CTmj  - Tot)* 


The  tube  metal  temperatures  were  obtained  from  an 
energy  balance/ 
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l**t  Eat  - E«*s 

for  the  screen,  and  for  the  main  bank 


2 . Desuperheater 

The  Dittus-Boelter  correlation  [12]  was  used  to 
determine  the  heat  transfer  from  the  steam  to  the  tubes,  i.e.. 


Eo*. 


(,LN\Tb} 


The  log  mean  temperature  difference  was 


LMTD  = ~ 

r (t.^  -TW.n 

The  energy  given  up  by  the  steam  was  obtained  from 


tot*  Cot  Grot  (Tbt_Tb3) 


The  heat  transfer  to  the  water  in  the  drum  was 


£«t 


-Tiwi) 


1 


J 


\ 

i 
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The  metal  temperature  was  obtained  from  the  energy  balance 


Cm  * E-oa  — 

d x 


3 . Economizer 

The  Dittus-Boelter  correlation  gives  the  heat 
transfer  from  the  tubes  to  the  water  as 

Grvoi  L*\t.6 


“Tw  j 


\T^  -T«i  J 


J 


The  heat  absorbed  by  the  water  was 

E.\m1  m Qr  Cw  ui  4 - Twul  } 

4 . Superheater 

The  heat  transfer  to  the  steam  was  given  by 

Gftv  * Ctqx  -tvO 

and  the  heat  transferred  from  the  metal  in  the  tubes 
was 
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T 


I 


tr6l  LV'Tfc 


l.w\to 


Tax  — Tbi. 

\J 

'MX-T6X 


The  metal  temperature  was  found  from 


Eex-Ebi 

dt 

E.  PRESSURE  DROP  CALCULATIONS 

For  turbulent  tube  flow,  the  pressure  drop  is  given 
by 

Al5  » V L/d 


Now 


Uwv  » 


Jr_ 

ft 


so 


AP 


A1  fU 


Lumping  the  constants  we  get 


; 


i 


i 


i 
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A?  * c 


£L 
f- 

Therefore  the  pressure  drops  in  the  superheater  and 
desuperheater  were  given  by 

T\»i  » Ct  pD1  • 

JO** 

and 

'P.i-f.j.  c, 

The  average  densities  were  found  from  the  ideal  gas  law 
under  the  assumption  that  the  temperature  variation  can 
be  neglected,  i.e.,  density  is  a function  of  pressure  only. 
This  assumption  is  motivated  solely  by  the  need  for  compu- 
tational efficiency.  It  is  reasonable  for  the  desuperheater 
since  the  temperature  excursions  ch^re  are  small.  Although 
errors  are  larger  for  the  superheater  density  values,  they 
are  still  within  accuracy  overall  limits  for  the  model. 

The  flow  through  the  throttle  value  was  given  by 

* Cj  T&x  Ve*t 
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F.  EQUATIONS  OF  STATE 


The  following  equations,  obtained  from  Reference  [13] , 
were  used  to  determine  state  properties.  Within  the  range 
of  pressures  for  which  the  model  is  valid  (300-1500  psi) , 
the  results  were  within  a few  percent  of  tabular  values  in 
the  worst  case. 

* E-XP  L C U fW  “ V \ 

TUv  - ExP  [ .iiisi  v 4nu3  \ j 

ft  “ M-B  - + l \37.£-5’Ts*.t  j 

— e-  a \ 

- 5^/^  - o.i  , 

- ‘t.i'ie-g'Psit 

j 
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IV.  DETERMINATION  OF  CONSTANTS  AND  INITIAL  CONDITIONS 

A.  GENERAL 

s For  the  model  to  be  useful,  it  should  be  possible  to 

obtain  the  necessary  constants  and  initial  conditions  from 
the  data  supplied  in  the  boiler  Technical  Manual  augmented 
with  reasonable  engineering  estimates  and  tabular  data  from 
standard  handbooks. 

The  primary  items  of  interest  are  heat  transfer  coeffi- 
cients and  initial  conditions  such  as  steam  flow  rates,  inlet 
, and  outlet  temperatures,  circulation  rates,  etc.  While 

some  of  these  were  directly  obtainable  from  the  boiler 

Technical  Manual,  others  such  as  tube  wall  temperatures  > 

and  riser  exit  quality  were  not.  Therefore,  a program  was 
developed  to  permit  the  user  to  enter  known  data  elements 
from  the  Technical  Manual  and  other  sources,  and  to  obtain 
as  an  output  the  necessary  values  for  the  constants  and 
initial  conditions  to  ensure  static  balance  at  the  speci- 
fied operating  point. 

The  basic  premise  under  which  the  program  operates 
* is  that  the  model  equations  with  all  derivatives  set  to 

zero  form  a determinant  set  with  the  addition  of  several 
relations  which  apply  under  equilibrium  conditions.  These 
relations  are: 

a.  Riser  and  downcomer  mass  flow  rates  in 
each  circuit  are  equal. 
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b.  Total  steam  flow  from  the  boiler  equals 
the  incoming  feed  flow  rate. 

c.  Total  heat  transfer  to  the  incoming  feed 
equals  the  heat  transfer  to  the  screen  and  the 
main  bank,  and  the  heat  removed  from  the  steam 
passing  through  the  desuperheater. 

d.  The  net  mass  transfer  at  the  steam  liquid 
interface  in  the  drum  is  zero. 

B.  PROGRAM  EQUATIONS 

The  total  heat  input  to  the  boiler  was  given  by 

■ HfifcT  ^avystc-^l1)  * Volume  (.  ) 

If  the  sensible  heats  of  the  oil  and  air  are  lumped  with 
the  heating  value  of  the  fuel,  the  heating  value  is 

Ho  * / Ga 

The  heat  absorbed  by  the  furnace  screen  at  steady  state 
was  calculated  from 

6.P.1  » Po«.HAct  HtAT  KftsoajVi'w  I see -■ft1) 
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Since  the  flue  gas  temperature  and  the  screen  tube  wall 
temperatures  were  provided,  the  heat  transfer  coefficient 
for  radiant  heat  transfer  to  the  screen  could  be  determined 
from 

/ { (Jk0!L)i-  ( r~JU)i}  . 

This  result  holds  for  the  conditions  reported  in  the 
manufacturers  technical  manual.  Degredation  due  to  age  or 
fouling  must  be  accounted  for  through  the  use  of  appropriate 
correlations . 

With  the  enthalpy  of  the  feed  water  at  the  economizer 

J 

outlet  and  the  enthalpy  of  the  steam  leaving  the  dry  pipe, 
and  the  change  in  enthalpy  of  the  steam  flowing  through  the 
desuperheater,  the  heat  transfer  to  the  main  bank  was 
calculated  from 

Hfcs  » Gro*.  ( ^ oi -eai  . 


The  metal  temperature  of  the  main  bank  risers  were  obtained 
by  calculating  the  constant 


and  solving  the  heat  transfer  equation 
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~Tlfca  * ^ Ea.3  / £**4  *?oi  ^ ^ ^-*-Tdi 

The  specific  heat  of  the  flue  gas  passing  the  main  bank 
was 

Cfcj*  " / G-fi.  C Too  — Tiiv ) 

and  the  unknown  constant  in  the  gas  side  heat  transfer 
relation  was  determined  from 

« £*3  / (Tfco. 3 -T«»)  . 

In  a similar  manner,  the  specific  heat  of  the  furnace 
flue  gas  is 

Cat  - (Eo-Eat')/  Ga.Tai 

and  the  water  side  heat  transfer  coefficient  is 

* Eai.  / ?0^‘  (rM i -Toi  ) * . 

The  heat  transferred  in  the  desuperheater  was: 
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With  this  value  and  the  temperature  difference  the  specific 
heat  of  the  steam  was  computed  from 


* Edi-  / 6-OL  (Tbt 


The  heat  transfer  coefficient  was  calculated  from  the 
Dittus-Boelter  correlation  111] 


The  log  mean  temperature  difference  was 
LftVTb  * 

and  the  metal  temperature  was  determined  by 
“I'm*  * (Tea  - QTo*  ) / ( \ - c9  } 


where : 


f=xp  [ ^ 


With  the  metal  temperature,  the  water  side  heat  transfer 
constant  was  obtained  from 

Kail.  • / CT at  — Twi  ") 


43 


r 


! 


i 


Analysis  of  the  Economizer  is  more  complicated  since  it  is 
a multipass  heat  exchanger  with  finned  tubes.  Additionally, 
since  we  are  concerned  with  energy  storage  within  the  metal 
of  the  economizer,  the  internal  and  external  heat  transfer 
must  be  considered  separately. 

The  heat  transferred  was  given  by  the  change  in  enthalpy 
of  the  feed  water,  i.e.. 


1 1 - £>oi 


The  specific  heat  of  the  water  was 

As  in  the  case  of  the  desuperheater,  the  inside  heat 
transfer  coefficient  was  calculated  from  the  Dittus- 
Boelter  correlation  as 


From  this  the  LMTD  was  obtained  as 

Un-*«  f i Cut 


I 
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Finally  the  metal  temperature  was  determined  from 


T"“\  “ -p-ru,3)/0-fj) 

P * «*p  { l«to  /Ct„m  _-rw,>  \ 


The  outside  heat  transfer  constant  was 

- Ewt  / Gr.  **  CTftAi  -Tmm') 


where 


"Tr.^3  * *Tte 

* 

The  specific  heat  of  the  flue  gas  passing  the  economizer 
was 


C*3  » Ewi  / Gt  CTaq  -TY^O 


The  superheater  heat  transfer  was  obtained  from  the  change 
in  enthalpy  of  the  steam,  i.e., 

E»i-  io-i-C  HeO 

With  this  value  and  the  known  metal  and  gas  temperatures, 
the  specific  heat  of  steam,  the  specific  heat  of  the  flue 
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gas  passing  the  superheater,  and  the  inside  and  outside 
heat  transfer  constants  were  determined  as 


CdI  s / Oat  CTot'Toi) 

Cftla,*  Eoi/Ga  CT(U  -Taj) 

- E01  / L^ro 

- Cb^/6a‘u  CToai.  - Tmi.) 

I 


where: 


LtftTO  * 


Tw.  — 
lTKV-Ttt».  \ 


and 


Toai  * CTat  +Ta»)  /x.o 


The  superheater  and  desuperheater  pressure  drop  constants 
were  evaluated  from  the  known  equilibrium  pressures  as 


C]  * ( *)  fo**  / Q0i.  1 

* CPox  Pb4LX  / Cox1 
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Similarly,  the  constants  used  in  the  calculation  of  the 
average  steam  densities  were  obtained  from: 


C / ( f0i 

and 


c'°'  to**'***}  . 

Finally  the  steam  value  flow  constant  was  computed  from 

£3  * 6dm  / £ pDi  l/nr) 

C.  CIRCULATION  BALANCE 

Since  they  are  not  normally  provided  by  the  manufacturer, 

values  for  the  friction  factors,  and  bend  and  entrance 

/ 

loss  coefficients  must  be  assumed.  The  values  do  not  appear 
to  be  particularly  critical,  and  tables  found  in  engineering 
handbooks  may  be  used  to  determine  reasonable  values  for 
these  parameters. 

An  iterative  procedure  must  be  used  to  adjust  the  various 
flow  values  involved  to  achieve  a steady  state  balance. 

Since  the  riser  exit  quality  for  boilers  of  this  type 
generally  lies  in  the  range  of  .02  to  .04,  a starting  trial 
value  for  the  circulation  flow  rate  may  be  obtained  from: 
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where  XE  is  the  exit  quality,  Er  is  the  heat  transfer  to 
the  riser,  Gr  is  the  flow  rate,  and  HFG  is  the  heat  of 
vaporization.  With  this  trial  value,  the  following 
balancing  procedure  is  used: 

1.  Set  HW2  (Downcomer  Enthalpy)  = HW1 

2.  Calculate  HWG  = HW2  + ED2/G0 

3.  Calculate  the  riser  exit  quality: 

Xo  - *9.\  / ( +■  VAu>i  — Hum*  ^ Gr» 

4.  Determine  the  non-boiling  height: 

W***0.  » (.^wt  - VtaO  / ( Hwl  + X0  VH}  — Www) 

5.  Compute  the  average  riser  density: 

6.  Obtain  the  outlet  density: 

*1'°/CVP  *•  v+0 
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7.  Compute  the  friction  factor: 


^ l i-iM  -aio^  1 ^ 1 * ,<'°18.s8  1,1 1 

8.  Next  the  riser  loss  factor  is  computed: 

9.  Finally  we  may  calculate  G0': 

Go  m ^ (.&>■»  ^ / C Cioo  v 

where  C100  is  the  downcomer  loss  factor.  ^ 

10.  Compare  Gjz 5 and  Gtf'  if  equal  go  on,  otherwise 
update  G0: 

Go  - Cro'-Go  + ^ 

* 

and  return  to  step  2. 

11.  Assume  a trial  value  for  Gs  (first  pass  only)  . 

12.  Compute  Xs : 

X*  ■ C&oi  — Go  Xo  ) /6*s 
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13-19  Same  as  for  the  main  bank  risers.  Return  to 


12  if  reiteration  necessary. 


20. 


Obtain  HW2 

VAiOj,  — 


from: 

Gun 
G*  *•  G’o 


— Hu*i ) +■  Hu»i 


21.  If  HW2 ' = HW2  go  on,  otherwise  update  HW2 : 


H*oi 


Hmi  — Mwr ^ 
Z 


and  return  to  step  2. 


Upon  exit  from  the  iteration  loop,  equilibrium  values 
of  the  flow  parameters  will  have  been  established.  With 
these  and  the  previously  calculated  heat  transfer  parameters, 
the  program  generates  the  list  of  constants  and  integration 
initial  conditions  required  by  the  dynamic  model.  The 
program  listing  and  a sample  output  is  given  in  Appendix  B. 
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V.  RESULTS  AND  CONCLUSIONS 


A.  OPEN  LOOP  RESPONSES 

The  output  curves  for  10%  increases  and  decreases  in 
throttle  valve  opening  and  fuel  flow  are  given  in  Figures  6- 
33.  Ramp  inputs  of  approximately  1.5%  per  second  were  used 
to  simulate  flexibility  test  conditions.  Only  one  input 
parameter  was  varied  for  each  run  with  all  others  held 
constant.  The  results  were  similar  to  those  presented  in 
[3]  and  [4]  insofar  as  general  trends  were  concerned.  Since 
the  boilers  involved  were  of  widely  different  geometries 
and  sizes,  a quantitative  comparison  was  not  possible. 

Since  CSMP  uses  only  single  precision  arithmetic  packages 
and  integration  algorithms,  the  high  frequency  noise  in  the 
outputs  is  attributable  to  the  roundoff  error  inherent  in 
the  wide  variations  in  magnitudes  of  numbers  used  in  the 
model,  e.g.,  quality  (order  of  0.01)  and  riser  mass  flow 
rate  (order  of  1000).  In  common  with  the  results  of  [3] 
and  [4] , the  "swell  effect"  in  response  to  an  increase  in 
throttle  opening  is  negligible.  It  should  be  noted  that  the 
swell  exhibited  in  the  data  of  Whalley  [3]  is  more  apparent 
than  real,  inasmuch  as  the  scales  in  his  graphs  are  very 
magnified  and  the  rise  is  actually  only  0.01  inches.  The 
sudden  increase  in  riser  mass  flow  rate  shown  in  [3]  is 
also  not  present  in  the  current  model. 
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Since  "shrink"  and  "swell"  effects  are  of  some  concern 
to  controls  designers,  an  attempt  was  made  to  model  for 
such  effects.  Four  possible  causes  for  the  phenomenon 
were  considered. 

1.  Homogeneous  Nucleation  in  the  Steam  Drum  Liquid 
This  effect  is  caused  by  a difference  between  the 

saturation  temperature  corresponding  to  the  actual  drum 
pressure,  and  the  actual  liquid  temperature.  The  super- 
heat required  for  such  nucleation  is  quite  high  however  and 
nucleation  will  occur  at  the  heating  surfaces  at  much  lower 
superheats.  Therefore  this  was  considered  to  be  an  unlikely 
cause  for  the  effect. 

2 . Nucleate  Boiling  at  the  Drum  Wall-liquid  Interface 
The  driving  force  for  this  phenomenon  is  also  the 

temperature  difference  between  the  drum  surface  and  the 
saturation  temperature.  Since  the  temperature  difference 
will  be  slight  except  under  large  drum  pressure  excursions, 
the  impact  of  this  factor  should  be  slight.  The  Rohsenow 
correlation  [5]  may  be  solved  in  terms  of  a rate  of  vapor 
generation  to  yield 


Cj  l3  / ft(fa-fv) 

c*  ^ Pt ] V v <r 


It  was  assumed  that  the  vapor  formed  would  be  uniformly 
distributed  bubbles  which  would  rise  to  the  surface  with  a 
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velocity  given  by  [14] 

Ub  ' iV  b 

\ ■ 

which  is  a conservative  estimate  since  interaction  effects 
tend  to  increase  the  rise  velocity.  The  bubble  diameter 
is  given  by  [14] 

T)  ■ *241  p ■,S" 

4 

This  gives  as  .25  ft/sec.  The  average  liquid  depth  is 
15  inches,  so  the  transit  time  is  = 5 seconds.  The  effect 
was  modeled  crudely  by  assuming  the  entrained  bubble 
mass  was  5.0(q/hfg)  lbm.  The  impact  of  this  model  on  shrink 
and  swell  values  was  negligible,  since  no  large  pressure 
excursions  were  developed. 

3.  Expansion  and  Contraction  of  Varpor  Bubbles  Entrained 
in  the  Riser  Fluid 

An  attempt  was  made  to  model  this  effect  by  including 
the  riser  tubes  in  the  control  volume  used  as  the  basis  for 
liquid  level  determination.  The  equation  used  was 

' 

v.l  = + (iSr-Ofc*  + 

where  Ralx  and  Ra2x  are  the  average  riser  densities  computed 
using  the  specific  volume  for  the  saturated  vapor  corresponding 
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to  drum  pressure,  while  Ral  and  Ra2  were  computed  from 
the  liquid  temperature  saturation  conditions.  The 
effects  were  insignificant. 

4 . Transient  Deviation  of  Riser  and  Downcomer  Mass 
Flow  Rates 

The  small  swell  effect  obtained  by  Whalley  [3] 
can  be  attributed  to  the  temporary  increase  of  riser  mass 
flow  rate  over  downcomer  mass  flow  rate.  The  period  over 
which  this  imbalance  could  be  sustained  is  short  and  would 
fail  to  account  for  the  longer  duration  of  "shrink"  and 
"swell"  effects  observed  in  practice.  It  may,  however 
trigger  one  of  the  two  phase  flow  instabilities  common  to 
natural  circulation  loops  e.g.,  flow  oscillations  [10,13]. 

The  analysis  of  these  effects  was  not  exhaustive,  and  the 
modeling  attempts  were  crude.  In  the  case  of  the  first  two, 
however,  tne  contribution  to  swell  is  probably  not  signifi- 
cant. Modeling  difficulties  with  the  latter  two  effects 
precluded  any  specific  conclusions. 

B.  CONCLUSIONS 

As  with  all  simulations,  the  information  available  from 
the  model  is  extensive.  Notably,  in  this  study,  main  bank 
and  water-screen  circulation  values  can  be  separately 
evaluated.  Hence,  the  contribution  of  each  of  these  ele- 
ments to  overall  boiler  operational  response  and  efficiency 
can  be  determined. 

The  model  simulates  the  responses  of  a Naval  boiler 
quite  well.  Unfortunately,  open  loop  data  is  unavailable 


54 


for  the  actual  boiler  and  insufficient  information  was 
available  to  develop  a model  of  the  existing  LHA-1  control 
system  as  a basis  for  comparison  studies.  As  such,  however, 
the  model  forms  a basis  upon  which  future  studies  involving 
modern  control  system  designs  can  be  established.  Moreover, 
while  considerable  refinement  of  the  model  is  possible, 
the  advantages  of  greater  accuracy  must  always  be  Weighed 
against  increased  computation  costs.  In  this  regard,  the 
implementation  of  more  efficient  integration  algorithms 
would  be  particularly  beneficial. 
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FIGURE  1.  SIMPLIFIED  BOILER  SCHEMATIC 
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FIGURE  19.  AVERAGE  RISER 
SETTING 


FIGURE  20.  WATER  LEVEL  RESPONSE 
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FIGURE  22.  STEAM  FLOW  RATE  RESPONSE  TO  A 10%  INCREASE  IN  FUEL  FLi 
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RISER  EXIT  QUALITY  RESPONSE  TO  A 10%  INCREASE 
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APPENDIX  A 


THIS  PAGE  IS  BEST  QUALITY  FRAGliCAiJLL 
EROM  COPY  EURN I SHE!)  TO  DOC 


THIS  PROGRAM  CALCULATES  THE  CONSTANTS  A NO  INITIAL  CONDITIONS 
FCP  « US  NAVY  O-TYPE  BCILcF  olVGN  ThS  INPUT  DA T A rB  TAINFC  F Kr'w 
The  ECILEP  TECHNICAL  MANUAL,  STEAM  T/sei:S,  ANC  STANDARD  ENG  I N Sr  - - 
INC  HANDBOOKS.  THE  INPLI  DATA  APS  ENTERED  IN  NAME  u I S T FORMAT. 
NAMELIST  DESIGNATIONS  AND  INPUT  NrMCN ICe  FOLLOW 


INCCM: 

QFFNT 

TCI  STM 

SPHSYM 

CSHSTM 

PC«LM 

SHCT 

SHCP 

CSHCT 

OSHCF 

:CCN  IT 

ECCNCT 

A I 5 TMP 

DIL7MP 

A I 3 F L C 

TIL  FLO 

X C S A I P 

DRAFT 

TGiSSC 

T 0 A S5H 

TGASMB 

TGASSC 

TSCRN 

TSFHTR 

HR  F VOL 

r HA  S 


CPFRATING  FC INT (PERCENT  I ♦ 

TOTAL  STEAM  F LCW (L3/HR ) 

SUPERHEATED  STEAM  FLOW  (LB/HP.) 
OrSUFERHfAT^D  STEAM  FlCW  (LS/HP) 

CRUM  PFESSUCC  (PSIG) 

SUP-FHEAT-p  CUTLET  TEMPERATLRC  (DGG-F) 
SUPERHEATED  OUTLET  PRDSSLP'  (FSIGt 
CE$UF=RHEAT'r  CUTLET  T "iMP  E R A”UP  i (CEC-F) 
OESUFFFHf ATEP  CUTLET  FRESSLP'i  (PSIC) 
ECONCMI ZCR  FE  r C INLET  T-MFERATUP"  (CoG-c) 
ECONCMIZER  f=cr:  OUTLET  TEMPERATURE  (CCG-F> 
AIR  INLET  TEMPERATURE  TO  REGISTERS  (TEC-?) 
OIL  TEMPERATURE  <*T  BORN-P  INLET  (OEG-F) 
AIRFLOW  RAT'  ( LB/HP  I 
CIL  FLOW  PATC  (LB/HA) 

EXCCSS  AIR  (PERCENT)  * 


CRAFT  LOSS  (IN-H20) 
FLU-  CIS  Tt MPr  RATURE 
FLUE  GAS  Tcf/p  = siTj^  = 
FLUE  GAS  TEMPCFATURE 
flue  gas  temperature 


LEAVING 
L EAV  INC 
LEAVING 
LEAVING 

SCREEN  TUBS  WALL  TEMPTraTL3: 

SUPEFHEATHp  TueE  WALL  TEMPERATURE  ( c e g- F I 
HEAT  RFL  E A SE  (K3TU/HR/CU  FT  Fijp\jArc  VELUM;  ) 
FURNACE  HEAT  A ESOP3T ICN  (KETJ/HR/SC  FT  RADIANT 
AESCRCING  APEA) 


SCREEN  (C^G-F) 
SUP"PH2A7r.P  (CEG-F) 
MAIN  EANK  (OEG-F) 
ECONOM 1 Z TR  ( OcG-  F ) 

( DOG-  F ) 


H F A ‘ 


NCTE:  "PERCEM™  VALLES  ARE  TO  BE  ICCOFD  AS  DECIMALS. 


SSCRcEN 

0TL3SC 
L A V SC 
NTLBSC 
RH  £S SC 
MASSSC 


AVERAGE  INSIDE  01  AMS  TER  OF  SCREEN  TLBES  (IN) 

AVERAGE  LENGTH  CF  SCREEN  TL  BE  (FT) 

NUMBEP  OF  SCREEN  TUBES 

RADIANT  HEAT  AESORBING  AREA  OF  FURNACE  SCREEN  (SC  FT) 
TOTAL  WEIGHT  of  SCREEN  TLBES  (LB) 


S SPHTR 


A'  SASH 
MASSSH 


TOTAL  SUPERHEATER  HEAT  TRANSFER  ARFA  (SC  FT) 
TOTAL  WEIGHT  OF  SUPERHEATER  TUBES 


JMNEANK 


CTUBM8  AVERAGE  INSIDE  DIAMETER  OF  MAIN  BANK  TLB  : S (IN) 
LAVMB  AVERAGE  LENGTH  OF  MAIN  BANK  TUBE  (IN) 

NTL3MB  NUMBER  OG  MAIN  BANK  TUBES 

Mf  S SMB  TOTAL  WEIGHT  CG  MAIN  oANK  TLS^S  (LB) 

A R E AMR  TOTAL  Hc  A T TRANSFER  APEA  (SC  FT) 


J.ECCN 


DTL3EC 

NFASSE 

NTL3EC 

LTLBEC 

MASSFC 


iNSirr  ciAMrT'p  of  economizer  tub'S(in) 

NUMBER  0 F t i j f c PASSES 
NUMBER  OF  TUBES  PER  PASS 
AVERAGE  LENGTH  OF  ECONOMIZER  TUBE  <FT! 
TOTAL  WEIGHT  OF  ECONOMIZER  (LE) 


iCESPF 


89 





.i. 


onnonnnnnnnnonnOnnoniinnri nnonnOonnonnnnnonnnnnnnnnni  mnonnonnonnonnonon 


THIS  PAGE  IS  BEST  QUALITY  PRACTICABIiE 
COPY  FURNISHED  TO  DDC 


CI1.2CS 
NTlBCS 
N F A S S C 
ITLBDS 
AREACS 
MASSCS 


IMSITE  D I A Md^  F 0 F DESUPERHEATER  ryef  (IN) 

NUMbER  OF  DESUPERHEATER  TUBES  PER  PASS 
NUMBER  CF  CP5UFERH;AT ER  Tie'  °ASSSS 
ID  NOTH  OF  CESUPERhEaTER  TIES  (FT) 

TOTAL  HEAT  TRANSFER  AREA  CF  D E SUP = PHc ATT R (SC  FT) 
TCTAL  WEIGHT  CF  DESUPERHEATER  ASSEMetY  (13) 


tCFMCCP 

CTLECC  AVERAGE  DIAMETER  OF  CP.UM  CCWNCOM-R  TieES  (IN) 
LAVDO  AVERAGE  LENGTH  OF  CRUM  DOWNCOMER  TURF  (FT) 
NTlBCO  NUMBER  OF  CFUM  DOWNCOMER  TCe:S 


$ FCRCC  3 

CTl  St-C  AV5RAGF  INSICF  DIAMETER  CF  HEADER  CCWNCO  l.ER  (IN) 
LAVFO  AVERAGE  LENGTH  CF  HEADER  DOWNCOMER  (FT) 

NTL3HO  NUMBER  OF  HEADER  OCWNCOMERS 


JECILER 

CSTMDM  CIAMETER  GF  STEAM  CRUM  (INI  ( 

LSTMCM  LENGTH  Oc  <t--AM  DRUM  (FT) 

CWTRCM  CIAMETER  GF  W A T c R DRUM  (IN) 

LWTPDM  LENGTH  OF  WATd  DRUM  (FT) 

HNCRM  FLIGHT  OF  NORMAL  WATER  LEVEL  ABOVE  BENCH  MARK  (FT) 

HHC R HEIGHT  OF  HEADER  (SCREEN)  ABOVE  BENCH  mark  (FT) 

HWTRDM  HEIGHT  OF  WATER  DRUM  ABOVE  BENCH  MA  0 K (FT)  ; 

F L P VOL  FURNACE  VCLUM*  (CU  FT) 


NCTE: 

STFEFMC 

HSHCUT 
HD  SCUT 
HEC  IN 
HE C CUT 
( 1IKH2C 

( 1)PRH2C 

( 1 ) VSCH20 
( 2 )KSTM 
( 2 )PRSTM 
( 2 JVSCSTM 
R SHOUT 
PCSCUT 
RFUS 

NCTES : 


CHOICE  OF  BCN'CF  MARK  IS  ARBITRARY 


ECONOMIZER  CUTLET  BnTHALPY  (BTU/lBMl 
DcSUFrRHEATFR  CUTLET  -NT  H A L P Y (3TU/IEM) 

ENTHALPY  C R ECONOMIZER  FREE  I NL0 T (fiTU/LBM) 

ENTHALPY  OF  ECONOMIZER  FE  EC  OUTLET  (BTi/LBM) 

THERMAL  CONCUC  T I VI TY  CF  WATER  ( 3T U / F R . FT  . C E G- F ) 

PRANCTL  NUMBER  FJR  WATER 

KIN'MATIC  VISCOSITY  FJR  *ATER  ( L B V/ F i -SEC ) 

THERMAL  CONDUCTIVITY  FOP,  ST  BAM  ( B TU/ H P . FT . CE  G- F ) 

PkANCTL  NUMBER  FOR  STEAM 

KINEMATIC  VISCCSITY  FOR  STEAM  (LBM/FT-SEC) 

SUP i PHE ATE R CUTLET  DENSITY  (LBM/CU  FT) 

CSSUPFRHSATFR  CUTLET  DENSITY  (LBM/CL  FT) 

OENSITY  OF  FLUE  GAS  AT  TGASSC  (LBM/CU  FT) 

(1)  RVALUATED  AT  AVERAGE  ECONOMIZER  WATER  TEMFERAncr 

(2)  EVALUATED  AT  AVERAGE  C E SUPERHS A T S3  STEAM  TEMPERATURE 


SSES 

KSCSC  ROUGHNESS  RATION  (SANC  ECU  IVALEN’-  ) FCR  SCREEN  TUP'  S 

KSCMB  ROUGHNESS  PATIC  (SAND  EQUIVALENT)  CF  main  bank  TuE' 

KSCCC  ROUGHNESS  PATIO  (SANC  EQUIVALENT)  FCR  93UM  OCwNCCVR 

KSCHD  ROUGHNESS  R AT I C (SANO  EQUIVALENT)  FCR  F-ADER  CCWNCr*'-R 

ENT  SC  SCREEN  'rUPE  ENTRANCE  LOSS 

eENCSC  SCREEN  Tues  BENO  LOSS 

cNTMB  MAIN  BANK  FNTPANuE  LOSS 

BENCMB  MAIN  BANK  BEND  LOSS  FAdCF 

EM  C C CRUM  DOWNCOMER  ENTRANCE  LCSS  FACTOR 

BENCCD  DRUM  DCwNC  C M E R BEND  LCSS  R A CTOR 

•MHO  HEADSS  DOW NC pMCf<  ENTRANCE  LCSS  RAdOR 

bendhd  header  DCWNCDMER  dENC  LOSS  FACTOR 


implicit  PEAL (A-Z ) 
INTEGER  I 
INTEGER  J.K.M.N 
INTEGER  NPASSc 
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& 


c 


c 

c 


THIS  PAGE  is  BEST  WALIH  PBiCIlCABW 
KHOU  COFX  FUKCiiSHHJ  TO  DDC  - 


DIMENSION  T F20  ( 20 ) » CPA'S  (20) 

DEFINE  NAMELISTS 

NAVELIST/TNCCNO/CPPNT.TCTSTM,  SPHSTM,CSFSTM,PORUV,'HCT,'HCP»C'HCT  , 
ID  SHOP ,=CONI  T.ECCNOT, A IPTMP ,31 LF LO, A IR F LO, XCSA  ! R .CRAFT  .TGASSC . 
2TCASSR.TGAS  V3.TCASEC,  TSCPN,  TS?HTR.FRF\,CL,CHAS  .CILfvp 
NAVE  LI $T/SCREEN/0TUBSC ,LAVSC, NTU3SC»FFASSC,MASSSC 
NAVEL  !ST/SuhTR/APFASH,VAS$SH 

NAVELIST/M'ieANK/LTUBVe,  LAVMB.NTUflMB  , V A SS  Me , AP.  G A V E 
NAVEL  l$T/SC0’l/0TU8rC  , NPA S SE ,NTUB  :C  , LTL  E 3C  , HAS  Sc  C 
NAv^l  IST/C'EPH/CTUB^S  ,NTUBnS,NPASSD,nL30S,APr'ACS,MASSCS 
NAVGLl'T/DSVDCR/CTUBCC  , L A VDO , NT UB DC 
NAVEL  I ' T /F3RDCR  /OTUBHD.LA  VHD,  NTG6H0 

NAVEL  1ST/B3IL=R/CSTMCV,LSTMQM,DWT8DM ,LwTROM,FNCF V ,FFOR  .FWTRTV  , 
iFLRVCL 

NAVEL  ! ST/TH  = *MO/H  SHOUT,  HP 'OUT  ,hSC  IN  ,HEC3UT, 

IK F20, PR H2C, VS CH2C.KST  V,P°'TM,  VS C STM  ,R  'FOUT  .ROSCLT  , PFLLE 
NAvELIET/LOSScS/KSDSC ,KSCVB,K'DDD,KSOC  ,5NtSC  .EFNCSC.ENTVE, 
JBcNCVE.ENTOC, BENCDD, rNTHD,BLUOHC 
NAVELl$T/nC0Nl/TP10,TM3,TM23,TM30,TV4D,TMS0,GCPC,GSRC,C>CC,;;)<SC 
NAVEL  1ST/1NC0N2  /HM10.YW10  ,VRO  J.HW60  ,GW2,TW3,GL,0e,GD2,GC10,V:NT 
NAVEL  1ST/  INC0N3/'~RA  10,CPA20,GQ1  £,SA  1XC  ,RA2X3,GWFC  ,G«SC 
NAVELI ST/C3NST1/KR.KS ,KCR ,KDS, AS, LP ,CR ,LR1, AS,LS, CS,LS1,LCP , ACR 
NAVEL  1 ST/C 0NST2 /DOR, HNBR, COS, LOS, HN9S  ,A0S,MV1  , V V2  , VH3  , V V4  ,M V5  , FX S 
NAVEL IST/C3NST3/FXS , V R 1 , VR,VS , V NDP ,VCF , LO IS , LC  IS 

NAVELIST/C3NSTA/CP1,Cc1A,CR2,CR3,CC1,CC2,CW,CV,KC1,KD2,KP1,K‘-'2,KRE 
NAVELIST/C3NST5/KRA,Z1,Z2,Z3,ZA,25,Zc,27,Z3,H1a2,GC13,CZ2,C8,CS,C10 
NAVEL IST/C3NST6/C11  ,C7 ,CSL»KWl,Kw2,KW-,KWA»HC 

NA«SL I ST/0UTPLT/SR1 ,F F 2 , G P3 , E R4 , - 0 2 , G'j , G$ , X T , X S , P A 1 ,RA2,F0, RS, 
JRw.HWl  , FW2, FW2,FW4,HW6,HFG 


READ ( 5 , INC3N3 ) 
wR  ITE ( 6 , 1 NC CNQ ) 

REAO  <£  , SCREEN) 

WP  ITE ( 6 » 60 1 ) 

WRITE  l£, SCREEN) 

READ  ( 5 , SPHTR ) 

W R IT  F ( 6 , 601 ) 

WRITE  (£, SPHTR) 

8 1 AC  ( £ » MN8 A NK  ) 

WP  ITE  (6,601  ) 

WRITE  (6  , VNBANK) 

READ ( 5 ,ECGN ) 

WRITS (6,601  ) 

WRITE ( 6 , E C 0 N ) 

REAC( £,D2SPH) 

WRITE  (6,601 ) 

WRITEU.DESPH) 

Re  AC  ( £ , CRMDCR ) 

WP  ITE  (6 ,601  ) 
hPITE  (fc.ORMDCR) 

R EAD ( 5 , HDRDCR I 
WRITE  <6  ,601  ) 

WP ITE ( 6 ,HORCCR ) 

RE  AD ( £ , BO IL  FR  ) 

WP ITE (6 ,601  ) 

write(£,bc:ler) 

REAC(5, THERMO) 

W ft  ITE  (6  ,601  ) 

WRITE(6,TH£RV0) 

REAC  (£  ,LCS3ES  ) 

WRITE  (6 ,601 ) 

WP  ITE<6, LOSSES) 

CALCULATION  CF  PRELIMINARY  CONSTANTS 


P I *3 . 1 A 15S2  7 

G02.2 

cv«. 1 1 

CRA10=0.0 

DPA2C=0.0 

ZER0=C.C 

Vf n r®c  PPNT 


t 
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KR  *KSCEB 
KS*KSCSC 
KCR-K'COD 
KCS  = KSOEiD 
2 1=5NT  CC 
Z2*BENCCD 
Z 2 * EN  1 H B 
ZA=BENCEB 
Z5  *ENT  S C 
ZA=BENCSC 
Z7*;NTEC 
Z 6 *B  ENDED 
PC1  = PCRLMUA.7 

EKl=EXP(0.2fcA52*AL0G(  PCI)  -« A . A 6 7 C 6 ) 
HFG=922.15-C.AC516*P01+1.717E-0A*PC1**2-A.219E-C8*PD1**3 
EC  1 = EW  1 + HPG 

TC1=P>°(. 22151 *ALOG (PCI ) *A. 77122) 

G8*GIIFIC/3600. C 

GL=A  1RFLC/380C.C 

TR13=TGASSC 

T N 10  = 1£CPN 

TE2D=1SPETR 

TP2=TP10 

1R  2 = T C A S SE 

TR  A = T G ASEB 

IPS^TCASEC 

GC  1 = TCTSTM/ 3633.0 

GC1£*CC1*.0C01 

GDA=SFHSTE/2600.0 

TC2=SECT 

P C2=  $ EC  P+ 1 A . 7 

TC  2=DSECT 

PC2  = CSECP«-1A.7 

CC2*CSESTN/3607  .0 

7 E 2 = EC0M  T 

TUA=ECCNCT 

E B =HRF VOL *FURVC 1*1000. C/ 3 600.0 

EP1iFEA$*RhASSC* 1000 .C/26CQ.0 

LS1=LAVSC 

LS=NUBSC*LAVSC 

CS*CTLBSC*NTU8SC/12.C 

AS  = FHCTUBSC**2*NTUbSC/<  A .On  AA.O) 

VS=AS*LAVSC 

PE1=EASSSC 

EE2=EASSSE 

L R 1 = LA  VMB 

LR=NTLBy'd*L  AVM3 

DF  = OU  E*8*NTUBME/12  .0 

AR=PI ♦0TU6MB**2*NTUBKe/( A.0*1AA .0) 

VP  = ARnAVM8 

PN3=NASSPB 

HMA=MA  S SEC 

NN5=EASSDS 

LCR*LAVDO*NTUBDC 

LOS* l A VHD*NTUBHC 

CCR=CTUBCD*NTUBCD/12.C 

DCS*CTLBEID*N  TUBED/ 12. C 

ACP  = NTUBDD* (P  I* < CTURO C /2 A .0  I **2 J 

ACS  = MLEHC-»  (P  I*  ( CTUBEC/2A  .0  )**2  I 

VCR=(PI*DSTN:0M«*2/(  A.0*1AA)  )*lSTMDE 

VNCR*  <PI*CWTRCM**2/  ( A .C * 1 AA  I ) *L V. TROE 

Ek3*ESC  IN 

EVA*EECCLT 

GR  *GB  «0L 

RR*6  2.8-0.01781*T0H-1.132E-05*TCl**2-t  . 786E- 08 *TC  1 < *3 
VF*1/Rh 

VFG  = 52A.0/PC1-.  1 
PDl«l.C/( VFEVFG) 

RC  2 = P SECUT 
RC3*R  C SCUT 
T W 1 = TC  1 
GW2»GC  1 

ER1*PFLLE*PURV0L 
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THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
FROM  COPX  FURNISHED  TO  DDQ ^ 

CZ2*(  CSTMDM/12.C MUST  MOM*  (1.0/12.0) 

CALCULATE  LCWER  HEATING  VALUE  OF  FUEL  CIL 
H0*EB/GB 

CCMPLTE  STOICHICMETRIC  AIP/FUEL  RATIC 


C 

C 

C 

c 

£ 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 


CSL=  ( 1 .0-XCSA  IR  MGl/GB 

RACIATION  HEAT  TRANSFFP  COEFFICIENT 

KFl=EFl/(  ( (TR10-*A60.0  )/100.0)  < TM  ] C + 463 . 0 ) / 10C.C  ) * *4  I 

CCMPUTS  HEAT  TR ANSC c R TO  TH-  MAIN  TUBE  BANK 
EP2=CC1*(HD1-HWK)-GD2* (HSFOUT-HDSOUT  » - ER 1 

COMPUTE  MAIN  BANK  AVERAGE  METAL  TEMPERATURE 

K V* 4 = C l.C/1.  7822  C6MA  REAM  E 

TM23  = (ER3/(  KN4*FD1**(  4.0/3.0m**(l.C/2.0MTDl 

COMPUTE  MAIN  BANK  GAS  SICE  HEAT  TRANSFER  COEFFICIENT 
KP2  = EF;/( ( ( TP-  + TRA) / 2 . 0- TM30 ) *GR* * . A ) 

COMPUTE  SCREEN  EANK  WATER  SIDE  HEAT  TRANSFER  C CEFFICILNT 
K V*  2=E  F 1 / < (PCI**  (A.  0/3.0  )*(TM10-T01  1**3) 

HEAT  TRANSFER  TC  Thc  SUPERHEATER 
£R2=GC1*(HSHCUT-FD1 ) 


SUPERHEATER  GAS  SIDS  HEAT  TRANSFER  COEFFICIENT 
KR2  = ER2/<  <( Tk2  + TP3)/2.C-TM2DI*GP**.6  ) 

HEAT  TRANSFER  TC  THE  ECONOMIZER 
EF4-CC1*<  HwA-HWE  ) 

CCMPUTS  SPECIFIC  HEATS  OF  FLUc  GAS 
FURNACE 

CP1=  (EE-ER1  >/ ( (TR10-8C.0  MGR) 

SUPERHEATER 

CP  1A=  ER2/  (GR*(TP2-TP.3  ) ) 

CCMEP2/(GD1*(TC2-TDU  > 

MAIN  BANK 

• CR2  = EP3/(GR*(TR  2-TR4 ) ) 

EC  C N CM  I ZEP 

CR2=ER4/(G?*(  TP.  4-TR  5 ) ) 

SPECIFIC  HEAT  FOR  WATER 
CW=SF4/(GD1*(TW'.-TW3)  ( 

COMPUTE  AVERAGE  HEAT  TRANSFER  COEFFICIENT 
CECUCTU6EC/12.C 

AREAEC=NTUdEC*L  TUBEC*PI*CECMNPASSE 

KWl=AFSAEC*<KH2C/DcCl*»l . S )*.  02  3*  ( 4 . C / l P I * VSC.H  2 C*NTU8EC  ) ) ** . 8 * 
JPFH2C**  .4 

LMTCEC=cR4/(KWl*GDl**.6) 

XXXX11=SXP< (TW4-TW31/LMTCEC) 

T M AO  - ( T W3-T  W4*XXXX11 ) / < 1 . O-XXXX  1 1 ) 

KP4=ER4/(GR**.6*( (TR4+TR5)/2. 0-TM401) 

SPECIFIC  HEAT  OF  STEAM  AT  AVERAGE  DESUPERHEATER  TEMPERATURE 
EC2MGC2MHSHOUT-HDSCUT)  ) 

CC2=EC2/(Gn:*(SHOT-OSHCT)  I 

DESUPERHEAT ER  S T CAM  SICE  HEAT  TRANSFER  COEFFICIENT 
0CSHMCTU80S/12 .0 

HXFRCS  = (KSTM/0CSH1**1 .8) * .023*1  4.0*CC2/<  P MVSC  STM*NTUBC$) )** .8* 
tPPSTN«*.3 

COMPUTE  THE  UCG  MEAN  TFM  PEPATGRc  CIFFSRENCE  FCR  THE  CE S U f ERH E AT  £ R 
LMTCCS=EC2/<HXFRCS*APEACS) 

DETERMINE  THE  DESUPf RHEATEP  METAL  TEMPERATURE 
XXXXC2=SXP(  ( SHOT-PjHOT  J/UMTOOS) 

TM  5 0 = ( S HGT-XX XXC2  *DS  HCT ) / (1 .0-XXXX02 ) 

OESUPERHFATER  WATER  SICE  HEAT  TRANSFER  COEFFICIENT 


KW2  = EC2/(TM50-Twl ) 
CESUPERHEATCR  STEAM  SIC= 


heat  transfer  coefficient 


KC2  = £C2/<  G02**.B*LMTCCS) 

superheater  steam  sicc  heat  transfer  CCEFFICIENT 


UMTOSH  = (TC2-TCI 1/AUOG ( (TM20-TD1  ) / ( T M 20 -TD2 ) ) 


KC  ME3  2/I  COM*  . E *UMTOIF  ) 

C ALCUU  ATE  THE  THFOTTUE  VALVE  FuCW  CCNSTANT 


nnn  o non  o moo  «~»ooor'*  nnonnnn 


i 


C11=0C4/(P02*VENT) 

C SLPERfEATER  OUTLET  DENSITY 

C8  = ( FC1+PP?  1/IRCI  + PC2) 
C10=(FC2+PD3)/(RL2*RD3> 

R A V SF  * l PC 1»  R02 ) / 2 . 0 
RAVCS = (RC2+SU3 I /2.0 
C7MFC1-PD2)*PAVSH/GC1**2 
C9  = < FC2-PD3) *RAVCS/G 0 2**2 


CIRCLLATICN  SYSTEM  CALCULATIONS 


COMPUTE  LOSS  COEFFICIENTS  FOR  DCWNCCNERS 

FCR=  < 1. 0/(1 .74-2.0* ALCG 10  (KS00D)) 1**2 

C1CJ  = <FCR*LCR /DC3  +Z 1*2  2 + 1 .0)/ (2.0*ACR**2*RW) 

FCS= ( 1 .0/(1 .74-2.0*ALCG19 (KSUhC  » ))**2 
C1C1=(FCS*LDS/DCS*Z7*Z6+1.C) /<2.0*ACS**2*RWi 

ITERATIVE  PROCcCURE  FCP  BALANCING  CIRCULAT!CN  LCCFS 

ScT  INITIAL  TRIAL  VALUES 

GC--2CC.C 
G S * 200  .0 
h V>2=H  H-EC2/G0 

START  FIRST  I TE  RAT  I CN  LOOP 

1 M»£-Fh2*ED2 /GO 

XC  = < ER2-G0*(HKl-HW6)  ) / < G0*HFG ) 

FXF  = FNCRM-FlhTRCN 

HN6R=FXR*AMAX1  ( <hhl-HW6)  ,0.0)/(  (HW1*XC*FFG)-Ffc6) 

L Xfi  = I-  X R-HNBR 

R A 1 = ( (LXR/  ( X0*V  FG  I ) * A L CG  ( ( XO  *V  F G ) / V F 1 . 1 >+RW*FNEP  ) /HXR 
RC  = i.C/(VF*XD-‘VFG) 

FR=(1.0/(1.  74-2 .0*4  LOG  10  ( KSDMB)  ))**2*(1. 0*1 6078. 58*(X0/FCl)**.Sfcl 
. PSIMB*  (FP.*  LR/CR-*1.0)*FA1/  (2. 3 ♦AR**2*RC*  *2  ) *Z  3/(2. G*AR**2*RW)+ 
tl 4/(2.C*AP**2*RC) 

G0X  = SCRT ( (KW-RA 1 ) *HXR*G/ (C100*PSIMB  ) ) 

400  FCP“AT(1H0,G15.5) 

IF  CFSCK  VALUE  EQUALS  INITIAL  GUESS,  CONTINUE;  ELSC  REITEPATE 
IF  (A6S ( GO-GO X ( .LT  ..01  IGO  TO  2 
GC*( GCX-GC) /2.C+G0 
GC  TO  1 

START  SECONC  ITERATION  LCCP 

2 CONTINUE 
XS=(GC1-G0*XC) /G  S 
FXS=FNORM-HFCR 

FNBS=HX$*A'/AX1  ( (HU-HW2)  . 0 . 3 1 / ( ( HW 1 ♦ X S * FFG  ) -HW  2 ) 

L X S = F X S-FN3  S 

RA2=((LXS/(XS*VFG))*ALCG((XS*VFG)/VF+1.0)+RW*FNES)/(HXS) 

RS*1.C/(VF*XS*VPG) 

F S=  ( 1 .0/  ( 1 . 74-2  .C*ALSG10<  KSOSC  m**2*(I.0*16D78.5E*(XS/FCl)**.96I 
PS  ISC  = (FS*L S/OS -*1. 0 )*PA2/  12.0 *AS*»2*PS**2)  *Z5/(2. 0*4 S**2*RW  I + 

$ Z 6 / ( 2.C*AS**2'*PJ  ) 

GSX=SCRT((RU-RA2)*HXS*G/«C101+PSISC)  ) 

IF  CHECK  VAUUH  EQUAUS  INITIAL  GUESS,  CONTINU6;  ELSE  REITEPATE. 
IF(ABSIOS-GSX>.LT..D1)CO  TO  3 
0S«(CiX-GS)/2.O4GS 
GC  TC  2 

START  THIPO  ITERATION  LOOP 

3 HK2X  = (GW2/<  GS*GC)  ) * ( H V4-HV, ! ) *hw  1 

C IF  CHECK  VALUE  EQUALS  INITIAL  GUESS,  CONTINUE;  ELS«=  REITERATE. 

I F (ABS  (HW2-FW2X  ) .LT.  . 1 )G0  TO  4 
H W 2 = ( F K2X-HW2  )/2  . 0*HW2 
GC  TC  1 

4 CCNTINLE 
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CALCULATE  IMTIAL  CCNCITICNS 

oorc  = cc  /p  c 

GSc0  = G£/P  s 

MWlC=VCR*Rw/2  .0 

HR10=RW1*MW10 

CXC3=RWl+X0*HFG/2.O 

CX$0=RWl+XS*hFG/2.0 

VRCJ  = (VDR/2.0) *Rri 

HW60=RW6 

GC  10=GC1 

GWR0=C0 

G W S0  = G S 

RA1XC  = (RAl  *HXR-RW*HKBF ) /LXR 
RA2X0-  (RA2*RXS-RW*HNB$  ) / L X S 
L C 1R  = L A VDG 
LD  1 S = L A VHD 
WRIT? (6,60)  ) 

6CC  FCFMAT(lHl) 

WRITE(6,INCCM) 

«R  ITS  (6,601  ) 

6C1  FCRVAT(IHO) 

WKITc(6,INCON2) 

WF  ITE  (6,601  ) 

WRITE  (t  .INCCN3) 

WRITE<6,601  ) 

WRITE (fc.CCNSTl  ) 

WRITE (6 ,601) 

W R IT  E ( 6 , CON  ST  2 ) 

WP  ITE (6  ,601  ) 

W R ITE  (6 ,CQN  ST3 ) 

WRITE (6,601  I 
WRITE (6.CCNST4) 

WRITE(6,601  ) 

WRITE  (6 , CON  ST  5 ) 

WPITE(6,fc01 ) 

WRITE (6 .CONST 6 ) 

WR  ITE (fc  ,600  ) 

WRIT?(6, OUTPUT) 

WRITE  <7 , INCONI ) 

WRITE ( 7 ,1 NCCN2  ) 

WRITE ( 7 , I NCCN3 ) 

W R IT  E ( 7 .CONST  1 ) 

W R I Tc ( 7 ,CCN  ST  2 I 
WR  IT  E ( 7 , CON  ST  3 ) 

WRITE  1 7 , CCN  ST  4 ) 

WRITE ( 7 , C C N S T 5 ) 

WP ITF ( 7 , CON  ST  6 ) 

S TCP 
ENC 
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